Segregation by thermal diffusion of an intruder in a moderately dense granular fluid 
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I. INTRODUCTION 
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A solution of the inelastic Enskog equation that goes beyond the weak dissipation limit and applies 
for moderate densities is used to determine the thermal diffusion factor of an intruder immersed 
in a dense granular gas under gravity. This factor provides a segregation criterion that shows the 
0^ , transition between the Brazil-nut effect (BNE) and the reverse Brazil-nut effect (RBNE) by varying 

the parameters of the system (masses, sizes, density and coefficients of restitution). The form of the 
phase-diagrams for the BNE/RBNE transition depends sensitively on the value of gravity relative 
to the thermal gradient, so that it is possible to switch between both states for given values of the 
. parameters of the system. Two specific limits are considered with detail: (i) absence of gravity, 

and (ii) homogeneous temperature. In the latter case, after some approximations, our results are 
\ consistent with previous theoretical results derived from the Enskog equation. Our results also 

indicate that the influence of dissipation on thermal diffusion is more important in the absence of 
gravity than in the opposite limit. The present analysis extends previous theoretical results derived 
in the dilute limit case [V. Garzo, Europhys. Lett. 75, 521 (2006)] and is consistent with the 
findings of some recent experimental results. 

■ Shortened title: Segregation in dense granular gases 
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One of the most important phenomena occurring in granular flows containing more than one species (a polydisperse 
system) is the segregation and mixing of dissimilar grains. This phenomenon, in which a homogeneous mixture of 
different species becomes spatially nonuniform by sorting themselves in terms of their masses and/or sizes, is of central 
I . interest in the field of granular matter mainly due to is industrial importance (powder metallurgy, pharmaceutical pills, 
glass and paint industries,. . .). The resulting non-uniformity is usually an undesirable property although there are 
some applications in which one wants to force species segregation (e.g., the separation of mined ores). Unfortunately, 
in spite of its practical relevance, the physical mechanisms that govern mixing and separation processes are not 
well understood yet. This fact motivates the development of fundamental theories that predict accurately the bulk 
behavior of these systems in order to be able to control such processes. 

It is well known that when a binary mixture constituted by one large ball and a number of smaller ones is subjected 
to vertical shaking in a container, usually the large particle (intruder) tends to climb to the top of the sample against 
gravity. This phenomenon is known as the Brazil-nut effect (BNE) and is one of the most puzzling problems of 
granular materials research [l|, |2|, y, On the other hand, a series of experimental works [5|, [6| have also observed 
£NJ ' the reverse buoyancy effect, namely, under certain conditions the intruder can also sink to the bottom of the granular 
, bed. This effect is known as the reverse Brazil-nut effect (RBNE). Several mechanisms have been proposed to explain 
the transition BNE/RBNE, for example, percolation [![, arching [H, convection 0, 0, 0]> inertia Q, condensation 
@, and interstitial-fluid effects Among the different competing mechanisms, thermal diffusion becomes the most 

■ relevant one at large shaking amplitude where the sample of grains resembles a granular gas. In this regime, binary 
^ j collisions prevail and kinetic theory can be a quite useful tool to analyze granular systems. 

The thermal diffusion factor has been recently @ evaluated for a dilute granular binary mixture from a solution 
of the inelastic Boltzmann equation that applies for arbitrary degree of dissipation and takes into account the non- 
equipartition of energy. The results show that the relative position of the large particles (intruders of mass mo) 
with respect to the small particles (gas particles of mass m) is determined by the sign of the control parameter 
(rrioT — mTo)/moT, where To and T denote the temperatures of intruder and gas particles, respectively. While in a 
molecular or ordinary gas mixture this sign is fixed only by the mass ratio of the particles (since Tq = T), for a granular 
gas mixture it also depends on the temperature ratio because of the lack of equipartition. This segregation criterion 
compares well with molecular dynamics (MD) simulations carried out in the tracer or intruder limit case 10]. The 
objective here is to extend the analysis made in Ref. Q to higher densities by considering the revised Enskog theory. 
By extending the Boltzmann analysis to high densities comparisons with MD simulations become practical and allow 
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one to quantitatively test the use of a hydrodynamic description for segregation in granular vibrated mixtures. In 
addition, at higher densities, it is possible that other segregation mechanisms different from the one identified in the 
dilute limit become relevant at those densities. 

Previous theoretical attempts based on kinetic theory to describe segregation in dense granular mixtures have also 
been reported. Thus, Jenkins and Yoon [Til ] have developed a hydrodynamic theory for the segregation of elastic 
particles, finding a criterion for segregation relatively close to the numerical results obtained by Hong et al. 0). Given 
that the criterion obtained in Refs. [6[ and (Tlj| only applies for elastic particles, more recently Trujillo et al. fl3j 
have derived an evolution equation for the relative velocity of the intruder by using the kinetic theory proposed by 
Jenkins and Mancini that only applies for nearly elastic particles. Interestingly, they considered the influence 
of the non-equipartition of granular energy (which is a generic feature of granular mixtures) on segregation through 
constitutive relations for the partial pressures. However, the results reported by Trujillo et al. [l3j have been derived 
by neglecting the presence of temperature gradients in the bulk region so that, the segregation dynamics of intruders is 
only driven by the gravitational force. Therefore, it appears that a complete theoretical description for the dynamics 
of BNE/RBNE in dense gases is still lacking. 

As said before, the goal of this paper is to analyze the segregation by thermal diffusion of an intruder in a dense 
granular gas. The segregation criterion is obtained from a recent solution [lilflU of the inelastic Enskog equation that 
covers some of the aspects not taken into account in previous works for dense systems [ill [l3j and extends previous 
results derived for dilute binary mixtures 0, to higher densities. Specifically, (i) it takes into account the nonlinear 
dependence of the transport coefficients on dissipation so that the theory is expected to be applicable for a wide range 
of values of the coefficients of restitution, (ii) it considers the combined effect of gravity and thermal gradients on 
thermal diffusion and (iii) it applies for moderate densities. Consequently, the theory subsumes all previous analysis 
for both dense and dilute gases, which are recovered in the appropriate limits. In addition, the theoretical predictions 
are in qualitative agreement with some MD simulations [ToL [la , [l7| and are also consistent with previous experimental 
works 18]. A preliminary report of some of the results presented here has been given in Ref. [13 ]. 

The plan of the paper is as follows. First, the thermal diffusion factor A is defined and evaluated in Sect. HT1 by 
using a hydrodynamic description. This factor provides a convenient measure of the separation or segregation of 
species generated by a thermal gradient in a multicomponent system. Once A is expressed in terms of the pressure 
and the transport coefficients associated with the mass flux of impurities, these coefficients are explicitly determined 
in Sect, imi bv solving the Enskog-Lorentz kinetic equation by means of the Chapman-Enskog method. This allows us 
to get A as a function of the parameter space of the problem, namely, the mass and diameter ratios, the coefficients 
of restitution for collisions among gas-gas and intruder-gas particles and the solid volume fraction. In Sect. IIV1 the 
form of the phase-diagrams BNE/RBNE is widely investigated by varying the different parameters of the system. 
Moreover, a close comparison with the theoretical results [ill, EH derived for thermalized dense gases is also carried 
out, showing that even in this limit the segregation criterion derived in this paper is more general than the one 
previously obtained since it covers the complete range of the parameter space of the system. The paper is closed in 
Sect. |V] with a brief discussion of the results obtained in this paper. 

II. HYDRODYNAMIC DESCRIPTION FOR SEGREGATION BY THERMAL DIFFUSION 

We consider a moderately dense granular gas of inelastic hard disks (d = 2) or spheres (d = 3) of mass m and 
diameter a. The gas is in the presence of the gravitational field g = — <?e z , where g is a positive constant and e 2 
is the unit vector in the positive direction of the z axis. The particles collide with a constant coefficient of normal 
restitution a. In the hydrodynamic description, it is assumed that the state of the gas is characterized by the local 
number density n(r,t), flow velocity U(r, t), and temperature T(r,t). The time evolution of these fields is given by 
the balance hydrodynamic equations 

D t n + nV ■ U = , (1) 
AU+(mn)- 1 V-P = g, (2) 



D t T + A (V • q+ PijVjUi) = -£T , (3) 
dn 



where D t = dt + u • V is the material time derivative. In the above equations, P is the pressure tensor, q is the heat 
flux and £ is the cooling rate associated with the energy dissipation in collisions. The macroscopic balance equations 
(H}-© are not entirely expressed in terms of the hydrodynamic fields, and thus do not comprise a closed set of 
equations. To close these equations one has to express the cooling rate and the fluxes as functionals of the fields. 
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Such expressions are called "constitutive relations" and provide the link between the balance equations and a closed 
set of equations for the hydrodynamic fields alone. Such a closed set of equations defines hydrodynamics in its most 
general sense. 

Let us assume now that some impurities of mass mo and diameter Co > o~ are added to the gas. Given that 
the impurities are present in tracer concentration, the problem is formally equivalent to study an intruder in a dense 
granular gas. This will be the terminology used in this paper. Collisions among intruder-gas particles are also inelastic 
and are characterized by the coefficient of normal restitution «o- It is also assumed that the presence of the intruder 
does not perturb the state of the gas and so, the flow velocity and temperature for the binary mixture composed 
by the dense gas plus the intruder are the same as those for the gas alone. Since the intruder may freely exchange 
momentum and energy with the gas particles, only the number density n (r,t) of the intruder is conserved. This 
continuity equation is given by 

D t n + n V ■ U + ^-^ = , (4) 

m 

where jo is the mass flux for the intruder, relative to the local flow U. 

In this paper, we are interested in analyzing segregation by thermal diffusion of the intruder in a granular dense 
gas [20( | . Thermal diffusion is caused by the relative motion of the components of a mixture due to the presence of 
a thermal gradient. As a consequence of this motion, a steady state is reached in which the separating effect arising 
from the thermal diffusion is balanced by the remixing effect of ordinary diffusion 20] . From an experimental point of 
view, the amount of segregation parallel to the thermal gradient can be characterized by the thermal diffusion factor 
A. Phenomenologically, it is defined at the steady state in the absence of convection (zero flow velocity) through the 
relation 



. <91nT 9 /no 

A— S — = IT ln — 

oz oz \ n 



), (5) 



where gradients only along the vertical direction (z axis) have been assumed for simplicity. Let us assume that gravity 
and thermal gradient point in parallel directions (i.e., the bottom plate is hotter than the top plate, c^lnT < 0). 
Thus, when A > 0, the intruder tends to rise with respect to the fluid particles (i.e., d z \a{no/n) > 0) while when 
A < 0, the intruder falls with respect to the fluid particles (i.e., d z ln(n /n) < 0). The former situation is referred to 
as the Brazil-nut effect (BNE) while the latter is called the reverse Brazil-nut effect (RBNE). 

As said before, we consider an inhomogeneous non-convecting steady state with only gradients in the z direction. 
Since U = 0, then the mass flux jo vanishes in the steady state according to the balance equation (j4]). Moreover, 
clearly the pressure tensor is diagonal for this state and so, Py — pSij where p is the hydrostatic pressure. In this 
case, the momentum balance equation ^fy reduces to 

^- = -pg, 6 
oz 

where p = mn is the mass density of the gas particles. As will be shown later, the spatial dependence of the pressure p 
is through its dependence on the number density n and the temperature T. As a consequence, Eq. ([6]) can be written 
more explicitly as 

E. dtT +^.d s n=-pg, (7) 

where the partial derivative d n p will be computed once the equation of state for the gas is obtained. To close the 
problem of determining A one needs a constitutive equation for the mass flux. Symmetry considerations yield 

ml dn m m dn p T dT 
p oz p oz T Oz 

where Do is the so-called kinetic diffusion coefficient, D is the mutual diffusion coefficient and D T is the thermal 
diffusion coefficient. These transport coefficients measure the contribution of each independent gradient to the mass 
flux of intruders. The condition jo ;Z = (applied to {SJ) along with the momentum equation Q allows one to express 
the thermal diffusion factor A defined by Eq. ([5]) in terms of the thermal gradient d z T and gravity. Its expression is 

_ (3D T * - (p* + g*)(D* a + D*) 

(3D* ■ W 
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Here, (3 — p* + nd n p* , p* = p/nT, and we have introduced the reduced transport coefficients 

D-* = A=D*, D* = ^D , ^ = =A (10) 
where r is a collision frequency (to be chosen later). In addition, 

■9* = JWr\ < (11) 

is a dimcnsionless parameter measuring the gravity relative to the thermal gradient. This quantity measures the 
competition between these two mechanisms (g and d z T) on segregation. 

To get the dependence of the thermal diffusion factor on the parameters of the system, the explicit form of the 
transport coefficients and the equation of state is needed. This can be achieved by solving the inelastic Enskog 
equation by means of the Chapman-Enskog method. This will be carried out in the next Section. 



III. TRANSPORT COEFFICIENTS 



We adopt now a kinetic theory point of view and start from the Enskog kinetic equation for the system. Thus, all 
the macroscopic (or hydrodynamic) properties of interest of the system (dense gas plus intruder) are determined from 
the one-particle velocity distribution functions of the gas /(r, v, t) and the intruder /o(r, v, t). Since thermal diffusion 
is given in terms of the transport coefficients Do, D, and D T associated with the mass flux of the intruder, our goal 
here is to solve the corresponding Enskog-Lorentz equation for the intruder by applying the Chapman-Enskog method 
pll ] to first order in the spatial gradients. The Enskog equation neglects velocity correlations among particles which 
are about to collide, but it takes into account the dominant spatial corrections to the Boltzmann equation (which only 
applies for dilute gases) due to excluded volume effects. Although the first assumption (molecular chaos hypothesis) 
can be questionable at high densities, there is substantial evidence in the literature of the accuracy of the Enskog 
theory for densities outside the dilute limit (moderate densities) and values of dissipation beyond the quasielastic 
limit. As a matter of fact, this is the only available theory at present for making explicit calculations of the transport 
properties of moderately dense gases. 

In order to fluidize the system in most of the experiments energy is added to the gas by the bottom wall which 
vibrates in a given way. Due to this external injection of energy, the system reaches a steady state whose properties 
far away from the boundaries (bulk domain) are expected to be insensitive to the details of the driving forces. 
However, due to the mathematical complexities associated with the use of vibrating boundary conditions, here we will 
introduce a stochastic external force, coupling the velocity of each particle to a white noise (stochastic thermostat). 
This kind of forcing, which has been shown to be relevant for some two-dimensional experimental configurations with 
a rough vibrating piston [22| . has been used by many authors [23j in the past years to analyze different problems, 
such as segregation in granular binary mixtures [13l. l l 91 ]. Although the relationship of these external forces with real 
vibrating walls is not clear to date, some results [9j derived in driven steady states for the temperature ratio by 
using the stochastic driving method agree quite well with molecular dynamics simulations [l6| of shaken mixtures. 
This agreement suggests that this driving method can be seen as a plausible approximation for comparison with 
experiments. In addition, the advantage of such a driving mechanism is that it lends itself to theoretical progress. 
Under the above conditions, the Enskog-Lorentz equation for the intruder reads 

(I + v • v 4S£ + g • I) ^ v < = j ^ ®> ™> ^ 

where the collision operator Jo[v|/o(t), /(*)] is 

Jo [n, vi | /o (*),/(*)] = a 4 - 1 J dw 2 J daG(a ■ g)(a ■ g) 

x K~ 2 Xo (ri,ri - CT)/ (ri,v";i)/(ri - a,v 2 ';t) 

-Xo (ri,ri +CT)/ (ri,vi;i)/(ri + W,v 2 ;t)} . (13) 

Here, <r = era, a — (<tq + cr)/2, <r is a unit vector along the centers of the two colliding spheres, ao (0 < ao < 1) is 
the coefficient of restitution for intruder-gas collisions, and xo is the pair correlation function for intruder-gas pairs 
at contact. The precollisional velocities are given by 

v'/ = Vl ^— (1 + a^ 1 ) (a ■ g)?, 

mo + m x 
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v£ = v 2 + — ^ (1 + a^ 1 ) (<r ■ g)£. (14) 
mo + m 

Upon writing Eq. (I12p wc have assumed that the system is driven by means of a stochastic Langevin force representing 
Gaussian white noise [24| . This force is written as ^o = ™o£, where the covariance of the stochastic acceleration has 
been chosen to be the same for the gas particles and the intruder [25|, [26|. In the context of the Enskog equation (|12[) . 
this external force is represented by a Fokker-Planck collision operator of the form —^((T/m)d 2 /dv 2 . Note that the 
covariance of the external force has been taken to achieve a constant temperature in the homogeneous state. The 
generalization of the force to the inhomogeneous case is essentially a matter of choice and here, for simplicity, we have 
assumed that the stochastic force has the same form as in the homogeneous case except that now £ and T are in 
general functions of space and time. This simple choice has been widely used in ordinary gases to analyze nonlinear 
transport in shearing systems [27j ■ 

As said before, the main goal of this Section is to compute the mass flux jo to first order in the spatial gradients. 
It is defined as 



Jo = m a 



J dvV/ (r,v,t), (15) 



where V = v — U is the peculiar velocity. At a kinetic level, another interesting quantity for the intruder is its local 
temperature defined as 



T ^-^J)j ^ 2 ^v,t). (16) 



This quantity measures the mean kinetic energy of the intruder. As confirmed by computer simulations [IcL [T5, l28l] . 
experiments [29j and kinetic theory calculations [301 ] . the global temperature T and the temperature of the intruder 
To are in general different. 

In order to determine the mass flux, we solve the Enskog-Lorentz equation by means of the Chapman-Enskog (CE) 
expansion [2l| . This method assumes the existence of a normal solution in which all the space and time dependence of 
/o occurs through the hydrodynamic fields no, n, u and T. The CE procedure generates the normal solution explicitly 
by means of an expansion in gradients of the fields: 

/o = / (0) + e/ (1) +---, (17) 

where e is a formal parameter measuring the nonuniformity of the system. The application of the CE method to the 
Enskog equation for polydisperse granular mixtures has been carried out very recently (l4l . [l5| in the undriven case. 
The extension of these calculations to the driven case is straightforward. We only display here the final expressions 
for the transport coefficients with some technical details given in Appendix lAl As for elastic collisions, the coefficients 
D$, D, and D T are given in terms of the solutions of linear integral equations, which can be approximately solved by 
considering the leading terms in a Sonine polynomial expansion. Here, we have considered for simplicity the first Sonine 
approximation. In dimensionless form, the transport coefficients are defined by Eq. (fTQ|) with t = na d ~ 1 y^2T/m [3lj ] - 
Their explicit forms have been obtained in Appendix lAl with the result 

D* Q = 4, (18) 

V D 

D T * = -™ (p* -1)+ {1+Uj)d — X °V(1 + ao), (19) 



T,n 



2v* 1 + M T V d<j) 



Here, 7 = Tq/T is the temperature ratio, M = mo/m is the mass ratio, u> = a^jo is the size ratio, Xo^ is the pair 
correlation function for intruder-gas evaluated at zeroth order, 

^-'"V^- xf ,/^ (lt4 (21) 



dr(f) w l + MV m 
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FIG. 1: (color online) Plot of the ratio of the mean square velocities mTo/moT as a function of h = m(l + a,<y)/2(m + mo) for 
hard disks (d = 2), a = 0.95, (Jo /a — 1 in the case of a dilute gas (0 = 0). The solid lines are the theoretical predictions given 
by Eqs. ([22} and (|27[) and the symbols are MD simulation results obtained by Brey et al. [TJJ] for different values of the mass 
ratio: mo/m = 2 (triangles), 1 (squares), 0.75 (stars), and 0.5 (circles). 



and fiQ is the chemical potential of the intruder. When the granular gas is driven by means of a stochastic thermostat, 
the temperature ratio 7 is determined from the requirement [26L |32| 



7Co = MC 



(22) 



where £q is the (reduced) cooling rate associated with the partial temperature To. The expressions of the (reduced) 
cooling rates and £q in the local equilibrium approximation are given by [lij 



t ~ dT(d/2) 



(23) 



t dr(f) Vo-y l + MV m 



7 + M ' 
1 -2 7 (l + M) (1 + ao) 



(24) 



where is the pair correlation function for the granular gas evaluated at zeroth-order. Moreover, the (reduced) 
pressure p* is [H, [3^ | 

p* = l + 2 d - 2 X (O) 0(l + «), (25) 

where 

no- d (26) 



r 2 d - 1 dT(d/2) 
is the solid volume fraction. 

Before considering the dependence of the transport coefficients on dissipation, let us illustrate first the dependence 
of the temperature ratio 7 on the parameters of the problem. It must be noted that condition (|22p to determine 
the temperature ratio differs from the one derived in the undriven (free cooling) case [3fj| . where 7 is obtained by 
requiring the equality of the cooling rates, i.e., 

C* = Co- (27) 

Since the condition to determine the ratio Tq/T is different in the driven and undriven states, it is interesting to 
explore the similarities and differences between the temperature ratios in both situations. Note that, according to 
Eq. (pM)) . the dependence of £0 on 7 is through the ratio of mean square velocities 8 = mT /m T. In terms of the 
parameter 0, the conditions (|22p and (|27[) are cubic equations with a unique real, positive solution. In particular, the 
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FIG. 2: (color online) Plot of the temperature ratio To/T as a function of the size ratio ao/cr for inelastic hard spheres (d = 3) 
with <j> — 0.1, a — ao = 0.8 and when the intruder and the gas particles have the same mass density [mo/m = (ao/cr) 3 ]. The 
lines are the kinetic theory results obtained in the driven and undriven cases. 




behavior of the solution in the limit O — > for undriven homogeneous states has been analyzed by Santos and Dufty 
[35l[36|. where a change similar to a second order phase transition has been shown. In Fig. [TJ we plot as a function 
of the dimensionless quantity h = mil + ao)/2(m + m ) for hard disks (d = 2), a = 0.95, (Jq/(7 — 1 and in the case 
of a dilute gas (0 = and so, x' ) = x^ = !)• The theoretical predictions obtained from the conditions and 
([27)1 indicate that is a function only of the parameter h for given values of a and (To/ a [35j . Molecular dynamics 
simulation results obtained by Brey et al. [101 ] in an open vibrated granular gas for different values of the mass 
ratio have been also included. It is apparent that, for the range of values explored in Fig. [IJ the heating mechanism 
slightly affects the value of O since the theoretical curves obtained from the driven and undriven conditions yield quite 
identical results. Moreover, the agreement between theory (driven and undriven cases) and simulation data is very 
good, except perhaps for small values of h where the results obtained in the undriven case compare with simulation 
data better than those derived in the driven case. On the other hand, significant discrepancies between the results 
obtained with and without a thermostat for the temperature ratio T /T are observed in Fig. [51 where T /T is plotted 
versus the size ratio ao/cr for a moderately dense gas {<fi = 0.1). Here, the intruder and the gas particles are composed 
by spheres of the same material and therefore, the same mass density [i.e., mo/m = (ao/u) 3 ]. To evaluate the ratio 
To/T, one has to know the explicit forms of x^ an d Xo°' > - ^ n the case of hard spheres (d = 3), a good approximation 
for x^ is provided by the Carnahan-Starling form [37| 



1 - 



l 



y(0) = _ (28) 

x (1-0) 3 ' [ ' 

while the intruder-gas pair correlation function is given by [38[ 

X ° = 1+LU (1-0) 2 + (l+O,) 2 (1-0) 3 ' { ' 

It is apparent that the disagreement found in Fig. [5] illustrates the fact that the heating mechanisms affect in general 
nonequipartition even in the bulk of the system [391 ]. Although the lack of simulation data in Fig. [2] prevent us to assess 
the reliability of both theoretical predictions, it must be noted that a recent comparison Q between MD simulations 
of agitated binary mixtures [16| and kinetic theory results based on the driven condition (|22[) shows a good agreement 
in contrast with the predictions obtained from the undriven relation (|27p where significant discrepancies between the 
latter and simulation appear especially as the size ratio increases (see Fig. 1 of Ref. [9]). More simulation data are 
needed to make quantitative comparisons between theories based on homogeneously heated granular systems and 
boundary-driven problems in order to asses the reliability of the above theoretical predictions. 

Now, the dependence of the transport coefficients on dissipation is considered. According to Eq. (|20p. in order to 
get the explicit dependence of Dq on a one has to give the form of the chemical potential /iq- The expression for the 
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FIG. 3: (color online) Plot of the reduced mutual diffusion coefficient D*(a)/D*(l) as a function of the (common) coefficient 
of restitution a = ao for inelastic hard spheres (d = 3) with mo/m = cra/cr — 2 and three values of the solid volume fraction: 
cj> — (solid line), (j> = 0.2 (dashed line) and cj> = 0.4 (dotted line). 
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FIG. 4: (color online) Plot of the reduced kinetic diffusion coefficient Do(a)/D '(l) as a function of the (common) coefficient 
of restitution a = ao for inelastic hard spheres (d — 3) with mo/m = oo/o — 2 and three values of the solid volume fraction: 
cj> — (solid line), <f> = 0.2 (dashed line) and cj> = 0.4 (dotted line). 



chemical potential of the intruder consistent with the approximation ()29|1 is [4C 
/'o 



T = C 3 +mn -ln(l -</>) + 3ujy-- 



3uj z 



21n(l - 
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FIG. 5: (color online) Plot of the reduced thermal diffusion coefficient D T * (a)/ D T * (1) as a function of the (common) coefficient 
of restitution a = ao for inelastic hard spheres (d = 3) with mo/m = co/cr = 2 and three values of the solid volume fraction: 
4> — (solid line), 4> — 0.2 (dashed line) and (j> — 0.4 (dotted line). 
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FIG. 6: (color online) Plot of the thermal diffusion factor A versus the (reduced) gravity \g*\ for mo/m — ao/u = 2, (f> = 0.2 
and three values of the (common) coefficient of restitution a — ao. The left panel is for hard disks (d — 2) while the right panel 
is for hard spheres (d = 3). 



where C3 is a constant. Figures H] and [5] show the reduced coefficients D*(a)/D*(l), Dq(o)/Dq(1), and 
D T *(a)/D T *(l), respectively, versus the (common) coefficient of restitution a — ao for inelastic hard spheres with 
mo/m = ao /a = 2 and several values of the solid volume fraction : </> = (dilute gas), <f> = 0.2 (moderately dense 
gas) and tj> — 0.4 (dense gas). All the transport coefficients have been reduced with respect to their elastic values 
D*(l), -Dq(I) an d D T *(1). We observe that in general the deviation from the functional form for elastic collisions is 
significant, especially in the case of the kinetic diffusion coefficient Dq- This deviation becomes more important as 
the density of the gas decreases. Although not shown in the above figures, a comparison with the results obtained in 
the undriven case 14 If shows again quantitative differences between the results derived with and without a thermostat 
and so, as expected [39j,|42|, the latter does not play a neutral role in the mass transport. 



IV. THERMAL DIFFUSION FACTOR. PHASE DIAGRAMS FOR THE BNE/RBNE TRANSITION 

Once the form of the transport coefficients is known, the thermal diffusion factor A can be explicitly obtained when 
one substitutes Eqs. (pi])-® for D*, D T * and D*, respectively, and Eq. (H) for p* into Eq. ©. This gives the 
dependence of A in terms of the parameter space of the problem. This space is sixfold: the dimensionless gravity g* , 
the mass ratio mo/m, the size ratio cro/c; the coefficients of restitution a and ao and the the solid volume fraction <p. 
According to Eq. ([9]), A is a linear function of gravity g*, as is illustrated in Fig. [6] where thermal diffusion is plotted 
as a function of \g*\ for hard disks (d = 2) and spheres (d — 3). We have considered the system mo/m = cro/a = 2, 
(j) = 0.2 and several values of the (common) coefficient of restitution a = ao- In the three-dimensional case, x*- - 1 , xb^ 
and no are given by Eqs. (f2"51) - (f3T))) . respectively, while for the two-dimensional case, and Xo^ are approximately 
given by [l2| 

* m ^ ,31) 



(0) 

Xo = 



1 



1 -<)> 81+w (1-0) 2 ' 
For hard disks, the form of the chemical potential consistent with the above approximations is [1^ 



Mo 
T 



= C2 + In no — ln(l — 



1 



9(f> 



+ ln(I - <j>) 



1 



Mf - ln(1 - +\ 



(32) 



(33) 



where C2 is a constant, ft is apparent for the case analyzed here that the RBNE is dominant for hard disks in all 
the range of values of gravity (except for a — 0.5 where a change to BNE is expected for values of \g*\ larger than 2) 
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while the RBNE is only dominant at small \g*\ for hard spheres. Thus, as already noted in Ref. [lj|, for given values 
of ?7io/m, <Jo/a, a, ao and <j> there is a critical value \g*\ such that a transition BNE<^ RBNE (or RBNE<=> BNE) is 
observed for \g*\ > \g*\. We see that the value of \g*\ increases with dissipation. 

The condition A = provides the segregation criterion for the transition BNE/RBNE. Given that the parameter 
(3 = p* + 4>dj,p* and the kinetic diffusion coefficient Dq = j/vj-, are both positive, then according to Eq. ([9]), the line 
delineating the regimes between BNE and RBNE is obtained from the relation 



(3D T * = (p* + g*){D^ + D*). 

This condition can be written more explicitly when one takes into account the explicit forms ([1 
coefficients. The result is 



(34) 

- ([20"]) of the transport 



g*( 7 -M) + <t> 



d-Mg*M 



M 



1 



(l+uj) d M (0) 



l + M' 



M + -7 

(p* + g*)^pA-[3 



= 0, 



where 



(1 



\d<j> 



(35) 



(36) 



7> 



Equation (|3^|) contains all the information necessary to describe the segregation due to thermal (or Soret) diffusion 
of an intruder in a moderately dense granular fluid. The first term on the left hand side measures essentially the 
influence of the non-equipartition of the granular energy on segregation. This term vanishes in the absence of the 
gravitational force. The second and third terms are proportional to the solid volume fraction <fi and so, they account 
for the effects of density on thermal diffusion. These latter two terms vanish in the dilute limit case (<j> — > 0). The 
influence of each term on the segregation criterion (|35[) depends on the specific values of dissipation (which is for 
instance the main responsible for the energy non-equipartition), density, mechanical parameters of the system and/or 
(reduced) gravity. 

Before exploring the dependence of the parameter space on the form of the phase diagrams BNE/RBNE, it is 
instructive to consider some special limit situations. Thus, when the intruder and particles of the gas are mechanically 
equivalent (too = to, ctq — c> and cto = ot), then as expected the two species do not segregate. This is consistent with 



Eq. [or Eq. ([53])] since in this limit case D T * = 



other hand, in the case of a dilute gas 



0), then Eq. 



D* = and so, 
5| yields 



9 *( 7 -M)=0. 



A = for any value of a and <f». On the 



(37) 



In the absence of gravity, Eq. (|3T[) holds trivially so that, the intruder does not segregate in a dilute gas when g* = 0. 
This result is due to the failure of the first Sonine approximation to accurately describe this special situation since 
segregation would appear if one would retain higher order terms in the Sonine polynomial expansion (20l . [4l| . On the 
other hand, when \g* \ ^ 0, the solution to Eq. (|3T|) is 



T 



mo 

TO 



(38) 



This result agrees with recent results derived from the Boltzmann equation [T(|. Note that, due to the lack of 
energy equipartition, the criterion (|38| is rather complicated since it involves all the parameter space. The segregation 
criterion (|38[) compares well with MD simulation results for the case of the steady state of an open vibrated granular 
system in the absence of macroscopic fluxes [Toj . 

According to Eq. (135|) . segregation is driven and sustained by both gravity and temperature gradients. The combined 
effect of both g and d z T on thermal diffusion is through the dimensionless gravity g* < defined by Eq. (fTTj) . This 
parameter measures the competion between both mechanisms on segregation. To separate the influence of each one 
of the terms appearing in (I35|) on segregation, we consider now some specific cases. 



A. Absence of gravity (|<?*| — > 0) 



Lets us focus first on a system in which effects of the temperature gradient are assumed to dominate over gravity so 
that the latter can be neglected. This limit is usually considered in MD simulations (see, for example, the simulations 
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FIG. 7: (color online) Phase diagram for BNE/RBNE for <j) — 0.25 in the absence of gravity (g* = 0) for two values of the 
(common) coefficient of restitution a = qq. Points above the curve correspond to A > (BNE) while points below the curve 
correspond to A < (RBNE). The dashed line is the result obtained for a = 0.7 assuming energy equipartition (To = T). The 
left panel is for hard disks (d — 2) while the right panel is for hard spheres (d = 3). 
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FIG. 8: (color online) Phase diagram for BNE/RBNE for inelastic hard spheres for a = ao = 0.8, \g* 
values of the solid volume fraction d>. 



and three different 



carried out by Galvin et al. [17|). In this case, \g*\ — > and Eq. reduces to 



dp* (l+m) d M (0) 



' d<t> 



P* 1 



j + M 
M 



(39) 



Of course, this equation is trivially satisfied in the case of a dilute gas (0 = 0). Beyond the dilute limit, the influence 
of each term in (|39p is still intricate. As an illustration, Figure [7] shows the phase diagram in the {mg/m, ao/cr}-plane 
at a total volume fraction of <fi — 0.25 (moderately dense gas) and two different values of the (common) coefficient 
of restitution ao = a. It is apparent that, in the absence of gravity, the main effect of dissipation is to reduce the 
size of the BNE. This effect is more significant in the case of hard spheres than in the case of disks. We observe that 
in general the RBNE is dominant for both small mass ratio and/or large size ratio. In order to assess the impact of 
the non-equipartition of granular energy on segregation, we have also plotted the corresponding phase diagram for 
a = 0.7 but assuming that Tq — T. The comparison between both curves clearly shows the significant influence of the 
temperature differences on thermal diffusion in the absence of gravity. This is consistent with the recent MD-fmdings 
of Galvin et al. [ijj where they showed that non-equipartition driving forces for segregation are comparable to other 
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FIG. 9: (color online) Phase diagram for BNE/RBNE for a two-dimensional system in the absence of thermal gradient (\g*\ — > 
oo) at <j> = 0.4 and two values of the (common) coefficient of restitution a = Qo. The dashed and dashed-dotted lines refer to 
the results obtained by Jenkins and Yoon [Tl[ for elastic gases (a = 1) and by Trujillo et al. [13| for a — 0.8, respectively. 



driving forces for systems displaying moderate level of non-equipartition. Figure [H] illustrates the influence of the 
volume fraction on the phase-diagram for a three dimensional system at moderate level of dissipation (a — 0.8). It is 
apparent that the role played by the density is quite important since the range of size and mass ratios for which the 
RBNE exists increases with decreasing <p. 



B. Thermalized systems (d z T — ► or \g*\ — + oo) 

In this Subsection the opposite limit is considered, namely, a system in which the global temperature of the 
bed does not vary with height (d z T — > 0). In this case, the segregation dynamics of the intruder is only driven 
by the gravitational force. This is a quite interesting limit since this situation (gravity dominates the temperature 
gradient) can be achieved in the shaken or sheared systems employed in numerical simulations and physical experiments 
[a EH H, Ell . Under these conditions (\g*\ — ► oo), the criterion ([55)) can be written as 

1 + 2 d_2 x ( W<?f>(l + a) [l + 4>d<f, ln(^x(°))] T m ' 1 ' 

As said in the Introduction, previous theoretical attempts to describe this particular situation have been made 
independently by Jenkins and Yoon [llj for elastic systems and by Trujillo et al. [13| for inelastic systems. Both 
descriptions are based on a kinetic theory [r2| that is restricted to the quasielastic limit (a — > 1), although Trujillo et 
al. [1 31 ] take into account the effect of non-equipartition of energy on segregation. Their segregation criterion differs 
from Eq. (I40[) and is given by [l3| 

-, . (l+w) d (0)- — 

1 + L T^«VTo_mo (41) 



l + 2 d - 1 X (0) (f> T m ' 

which is consistent with the one derived by Jenkins and Yoon (llj | when a = ao = 1. The discrepancies between 
Eqs. (|40p and (|4"Tj) can be attributed to the simplicity of the kinetic theory used for deriving the latter criterion. In 
particular, it is easy to see that Eq. ([40]) reduces to Eq. (|4"Tj) when one (i) neglects the dependence on inelasticity 
and assumes equipartition in certain terms, (ii) takes the approximation A = 1 (which only applies for a dilute gas of 
mechanically equivalent particles), and (iii) neglects high density corrections (last term in the denominator of l|40p). 
Therefore, in relation to the above previous results, we can conclude that the criterion (|4T))) is much more general 
than the one derived by Trujillo et al. [H| in the limit \g*\ — -» oo since our results cover the complete range of the 
parameter space of the problem. 

A typical phase-diagram for thermalized systems delineating the regimes between BNE and RBNE is plotted in Fig. 
[5]for the two dimensional case. (The qualitative features of the corresponding phase-diagram for the three-dimensional 
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FIG. 10: (color online) Phase diagram for BNE/RBNE for inelastic hard spheres for a = ao = 0.7, \g*\ = 1 and three different 
values of the solid volume fraction (f>. 

case are similar). Comparison between the left panel of Fig. [7] (hard disks) with Fig. [9] clearly shows that the presence 
of gravity changes dramatically the form of the phase-diagram. In p articular, the main effect of inelasticity is to reduce 
the size of RBNE region, which is consistent with experiments [la ]. Moreover, we also observe that the RBNE regime 
appears essentially now for both large mass ratio and/or small diameter ratio. On the other hand, the predictions of 
Trujillo et al. (l3j for a = 0.8 disagree with our results even at a qualitative level since they find that the mass ratio 
is a two- valued function of the size ratio in the phase diagram. In fact, according to the results of Trujillo et al. [l3T ]. 
the effect of dissipation is to introduce a threshold size ratio above which there is no RBNE. We also observe that 
our results differ from those obtained by Jenkins and Yoon (TJ for elastic gases, especially for large size ratios. Our 
results also indicate (not shown in Fig[5]) that non-equipartition has a weaker influence on segregation for thermalized 
systems than in the opposite limit (\g*\ = 0). This behavior qualitatively agrees with the experiments carried out 
by Schroter et al. 16] for vibrated mixtures as well as with some recent theoretical results of Yoon and Jenkins [46[ 
since both works find that segregation (when is only driven by gravity) is not significantly influenced by the difference 
between the temperatures of the two species. 



Finally, we consider the effect of density for finite values of the reduced gravity |<7*|. Figure ITD1 shows a phase 
diagram when \g*\ = 1 (gravity comparable to the thermal gradient) for different values of the volume fraction. We 
have considered inlastic hard spheres (d — 3) with = a = 0.7. In contrast to Fig. we observe that the REBNE 
regime appears essentially now for both large mass ratio and/or small size ratio. Regarding the influence of density 
on the form of the phase-diagram, it is apparent that the regime of the RBNE decreases significantly with increasing 
(f>. Following Trujillo et al. [HI, in the fluidized regime the effect of shaking strength of vibration on the phase- 
diagram for BNE/RBNE can be tied to the effect of varying the solid volume fraction <f>. According to this argument, 
Fig. [TUl shows that the possibility of RBNE increases with increasing shaking strength (or decreasing density). The 
experimental findings of Breu et al. fl8T ] show similar trends with increasing shaking strength, which is consistent 
with our results. 



The problem of segregation by thermal diffusion of an intruder in a dense granular gas has been addressed in 
this paper. Thermal diffusion is the relevant segregation mechanism in agitated granular mixtures at large shaking 
amplitude. In this situation, the so-called thermal diffusion factor A characterizes the amount of segregation parallel 
to the thermal gradient [20j . Here, the factor A has been obtained in a nonconvecting steady state with gradients only 
along the vertical direction (parallel to gravity). Two complementary approaches have been followed to evaluate the 
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thermal diffusion. First, by using a hydrodynamic description A has been expressed in terms of the pressure of the 
granular gas and the transport coefficients associated with the mass flux of impurities. Then, the above quantities have 
been explicitly determined by solving the inelastic Enskog equation by means of the Chapman- Enskog method [2lj |. 
This allow us to determine A as a function of the mass and size ratios, the coefficients of restitution for collisions among 
gas-gas and gas-intruder particles, the solid volume fraction and the reduced gravity g* [a dimensionless parameter 
defined in Eq. (fTTj) ]. Once the explicit form of A is known, the condition A = provides the segregation criterion for 
the transition BNE<^>RBNE. This criterion is given by Eq. (|55|) and is the most relevant result of this paper. 

Although some previous theoretical efforts [111 . Il3j on the same topic for dense granular gases have been made, 
they have been based on a kinetic theory which is valid for nearly elastic particles and have considered situations 
where gravity dominates over the temperature gradient (and so, the effects of the latter on segregation have been 
neglected). The present study goes beyond the weak dissipation limit and takes into account the influence of both 
thermal gradient and gravity (through the reduced gravity g* — pg/nd z T < 0). In addition, previous results [9l. Hoj 
derived in the dilute regime limit are recovered at zero density (<j) — > 0) . 

To illustrate the form of the phase-diagrams BNE/RBNE in the mass and size ratio plane, two specific situations 
have been mainly studied: g* = (absence of gravity) and \g*\ — > 00 (homogenous temperature). The first case has 
been considered in recent MD simulations [17j | while the second case has been widely studied by using kinetic theory 
[HI [HI , computer simulations [l6[ and experiments [l8[ . Our results show that the influence of dissipation on thermal 
diffusion is more important when the thermal gradient dominates over gravity (g* — 0) than in the opposite limit 
(|<7*| — ► 00). This weak influence on dissipation in the latter case contrasts with the results of Trujillo et al. j 1 31 ] since 
they found the main effect of inelasticity is to introduce a threshold size ratio above which there is no RBNE (see Fig. 
H]) . We attribute this discrepancy with Ref. [l3[ to the use of some (uncontrolled) approximations in the expressions 
of the partial pressures and the transport coefficients. Regarding the role played by the non-equipartition of granular 
energy (pseudo-thermal buoyancy force) in the segregation process, our results indicate (see Fig. [7]) that while the 
temperature differences has an important influence on thermal diffusion in the absence of gravity, it has a weaker 
effect on segregation when gravity dominates over thermal gradient. These conclusions agree qualitatively well with 
recent MD simulations [l7j and with some experiments carried out by Schroter et al. [l6j in vibrated mixtures. 

Although the theory reported in this paper is consistent with previous numerical and experimental results, a more 
quantitative comparison with the latter would be desirable. As a first test, kinetic theory predictions in the Boltzmann 
limit ((f) — > 0) 9| compare well with MD simulations of agitated dilute mixtures [lOfl . Given that the results derived 
here extend the description made in Ref. [9j to moderate densities, it can be reasonably expected that such a good 
agreement is also kept at finite densities. In this context, it is hoped that this paper stimulates the performance of 
such simulations. Another important limitation of the present study is its restriction to the tracer or intruder limit. 
This precludes the possibility of analyzing the influence of composition on the thermal diffusion factor. The extension 
of the results derived here to finite mole fraction is an interesting open problem. Moreover, this study will allow us 
to compare the theory with the results reported by Schroter et al. [16| in agitated mixtures constituted by particles 
of the same density and equal total volumes of large and small particles. To the best of my knowledge, this is one of 
the few experiments in which thermal diffusion has been isolated from the remaining segregation mechanisms. 
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APPENDIX A: CHAPMAN-ENSKOG SOLUTION IN THE DRIVEN CASE 



In this Appendix we solve the Enskog-Lorentz equation (|12[) to first order in the spatial gradients by means of 
the Chapman- Enskog method [2l[. From this solution, we determine then the transport coefficients Do, D and D T 
associated with the mass flux. 



The zeroth-order distribution function /q obeys the Enskog equation 



C (0) t a 2 (o) 

2m 8V 2Jo 



where 



d V2 / dae(a ■ g )(? • g ) k 2 /j c) K)/ (0) K)/ (0) (v 3 ) 



(Al) 



(A2) 
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Upon writing Eq. (|A1[) use has been made of the fact that d^T = at this order in the driven case. Moreover, 
refers to the cooling rate evaluated by using the zeroth-order velocity distribution function /(°). Since the latter is 
not exactly known, one has to expand /^°^ in Sonine polynomials. A good approximation for it (at least for not very 
strong inelasticity) is given by the Gaussian distribution 

/ (0) W-/M(^) = n(^) d/2 e-^ 2 / 2 T. (A3) 



Employing it, one gets the expression (|23[) for the reduced cooling rate C* = C^°V T where r = na d 1 y / 2T/m. The 
lowest order in the gradients of the cooling rate £q ^ for the temperature of the intruder T can be also obtained by 



approximating fjf ^ by the Gaussian distribution 



C(V) - foMV) = n ^y 2 e -^ftT (A4) 

With this approximation one gets expression (|24|) for Q = Q /r. 

The determination of the first order distribution follows similar mathematical steps as those made in the 
undriven case for polydisperse systems [3]. Here, we only display some partial results in the tracer limit (no/n — > 0). 
The distribution /q can be written as 

-> A)(V) • VlnT + B (V) • Vlnn +C (V) • Vlnn 

+Vo,ij (V) i (diUj + djUi - ~<%V ■ u) + £ (V) V • U, (A5) 

where the quantities «4o, Bo, Co, T^o,ij and £q are the solutions of a set of coupled linear integral equations. In this 
paper we are only interested in the first-order contribution to the mass flux j . It is defined as 

j^=m o ydvV/ (1) (V). (A6) 

Use of Eq. (|A5|) into Eq. (|A6j) and taking into account symmetry considerations, one gets the constitutive form for 
jo given by Eq. |[8} where 



D 1 



is the thermal diffusion coefficient 



is the kinetic diffusion coefficient and 



Do 



-!^JdvV-Ac(V) (A7) 



J dvV ■ B (V) (A8) 



rrionod 



D = -^JdvV-Co(V) (A9) 



is the mutual diffusion coefficient. According to Eqs. (|A7[) - (|A9[) . only the coefficients A®, Bo, and Co are involved 
in the evaluation of the n 
linear integral equations: 



in the evaluation of the mass transport jg of the intruder. These quantities are the solutions of the following set of 



- 4 0) [^J (0) ] = Ao + J§» [/ (0) , A], (A10) 



2 



2m dV 2 



B -4 0) [Bo,f {a) }=B , (All) 



C {0) T d 2 „ „„ _ l „ ,„ 



2m dV 2 



Co 4 v '[Co, r>] = Co + jn/r,c]. (A12) 
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The inhomogeneous terms Aq, Bq and Cq of the integral equations (j A10|) — (|A12|) are defined by 



Ao,(v) = iy i A.( V /r 



(0) 



P dVi 



(A13) 



Bo (V) = -V/a 



(0) 



(A14) 



^o,i{v) — -m - — — 77y\ ^rr 

9n dVi x w T v y T no 



(0) 



(A15) 



In Eqs. (|A13|) - (|A15|) . the pressure p is given by [301134} 

p = nT [1 + 2 d - 2 X <f>(l + a)] , (A16) 
a = aa/a is the size ratio, /xo is the chemical potential of the intruder and the operator ZCo,i[A^ is defined as 

tC 0<i [X] = o*^ J dv 2 J dZQ(B ■ g)(<r • g)a, [a^ 2 C (V'I)X(V^ + / (0) (Vi)X(V 2 )] . (A17) 

Note that, in contrast to what happens in the undriven case (l4| . here each one of the quantities »4o, 2?o, and Co 
obey closed integral equations. Moreover, upon writing Eqs. (jAlOj) — (|A12|1 . use has been made of the expression of the 
first-order distribution of the gas particles. Its form in the driven heated case has been derived in Ref. [13] and 
reads 



/ 



(i) 



A(V) ■ VlnT + C(V) • Vlnn 
1 



(V) 2 ( diUj + djUi - -<%V ■ U ] + £ (V) V • U. 



(A18) 



For practical purposes, the integral equations (|A10[) - ||A12I) must be approximately solved by using a Sonine poly- 
nomial expansion. In the lowest Sonine approximation, the quantities «4o, So and Cq are approximated by 



A)(V)^-/ ,m(V)^-VC t , 

^0^0 



(A19) 



Bo(V)^-/ ,m(V)-^VA), 
pT 



(A20) 



m 



n To 



(A21) 



where /o,m is given by Eq. (|A4[) . Consistently, A, and C must be also approximated in a similar way. However, 
both quantities vanish in the lowest Sonine approximation [331 ] - To get the transport coefficients D T , Dq and D, we 
substitute first *4q, Bo and Co by their Sonine approximations (j A19|) — (|A2L|) . respectively, and then, we multiply the 
integral equations (|A10[) - (|A12[) by nioV and integrate over the velocity. After some algebra, one gets 



vdD 1 



rap 



1 - 



m p 



2dp 



dV TOoVi/Co... 



w ■ ( v ' l0) ) 



(A22) 



n P T ° 



(A23) 



v D D 



xoMn dp (1 + uj) d ( dpo 



mo dn 



d X ( o 0) T 



) I dV ViKo 

T,n 



f(0) 



(A24) 
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where xq — ng/n, M = mo/m is the mass ratio and 

1 



v D = 



dn T 



dvmoV-4 0) [/o,MV,/( )] 



(A25) 



Note that Eqs. (|A22[) - (|A24[) have been obtained by neglecting some non-Gaussian contributions to the zeroth-order 
distribution /q°\ The collision integrals appearin g on the right hand side of Eqs. (|A22|) and (|A24[1 along with the 
collision frequency Vry have been evaluated in Ref. [151 ] in the case of a multicomponent mixture. In the tracer limit, 
one easily gets 



— / dV m Q VilC 0t i 



dV V J 



2 m y ' 1 + M 



- j dV m ViJC . 



f(0) 



1 7V1 _L \d l + M j. (0) 



(A26) 
(A27) 



VD = — — ?tcr 



dT(|) 



(0) 



m 1 + M 



7 + M 
M 



(1 + ao), 



(A28) 



where 7 = Tq/T is the temperature ratio. The final expressions for the reduced transport coefficients defined by (|10| 
can be easily obtained from Eqs. (|A22|> — (|A24ji and (|A26|) — (|A28|) . They are given by Eqs. (fl"8 l) -([20] ) with v* D = v d /t 
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